Molecular Modeling of the Adsorption of an Egg Yolk Protein on a Water-Oil Interface.

Egg yolk contains several molecular species with emulsifying properties, such as proteins and phospholipids. In particular, these molecules have both polar and nonpolar parts and thus can act as surfactants. One of the most surface-active proteins from egg yolk low-density lipoproteins is the so-called Apovitellenin-1. Experimental studies have been hindered by difficulties in isolating individual species from egg yolk lipoproteins. The purpose of this work was to assess the emulsifying properties of Apovitellenin-1 and any potential cooperative or competitive behavior in the presence of phospholipids. To do so, molecular simulations were carried out in a liquid-liquid interfacial system consisting of water and soybean oil, with varying concentrations of phospholipids and for different spatial configurations. To evaluate the conformational stability of the protein at the water-oil interface, the Gibbs free energy was computed from Metadynamics simulations as a function of the distance from the interface and of the radius of gyration. Moreover, a detailed analysis was also performed to determine which peptide residues were responsible for the protein adsorption at the oil-water interface as well as the lowering of the interfacial tension. Lastly, we combined the simulation results with a thermodynamic model to predict the interfacial tension behavior at increasing protein bulk concentration, which cannot be measured experimentally.

Mechanistic Investigation of tert-Butanol's Impact on Biopharmaceutical Formulations: When Experiments Meet Molecular Dynamics.

The use of tert-butyl alcohol for the lyophilization of pharmaceuticals has seen an uptick over the past years. Its advantages include increased solubility of hydrophobic drugs, enhanced product stability, shorter reconstitution time, and decreased processing time. While the mechanisms of protein stabilization exerted by cryo- and lyo-protectants are well known when water is the solvent of choice, little is known for organic solvents. This work investigates the interactions between two model proteins, namely, lactate dehydrogenase and myoglobin, and various excipients (mannitol, sucrose, 2-hydroxypropyl-ß-cyclodextrin and Tween 80) in the presence of tert-butyl alcohol. We thermally characterized mixtures of these components by differential scanning calorimetry and freeze-drying microscopy. We also spectroscopically evaluated the protein recovery after freezing and freeze-drying. We additionally performed molecular dynamics simulations to elucidate the interactions in ternary mixtures of the herein-investigated excipients, tert-butyl alcohol and the proteins. Both experiments and simulations revealed that tert-butyl alcohol had a detrimental impact on the recovery of the two investigated proteins, and no combination of excipients yielded a satisfactory recovery when the organic solvent was present within the formulation. Simulations suggested that the denaturing effect of tert-butyl alcohol was related to its propensity to accumulate in the proximity of the peptide surface, especially near positively charged residues.

Force Field Parameterization for the Description of the Interactions between Hydroxypropyl-ß-Cyclodextrin and Proteins.

Cyclodextrins are cyclic oligosaccharides, widely used as drug carriers, solubilizers, and excipients. Among cyclodextrins, the functionalized derivative known as hydroxypropyl-ß-cyclodextrin (HPßCD) offers several advantages due to its unique structural features. Its optimal use in pharmaceutical and medical applications would benefit from a molecular-level understanding of its behavior, as can be offered by molecular dynamics simulations. Here, we propose a set of parameters for all-atom simulations of HPßCD, based on the ADD force field for sugars developed in our group, and compare it to the original CHARMM36 description. Using Kirkwood-Buff integrals of binary HPßCD-water mixtures as target experimental data, we show that the ADD-based description results in a considerably improved prediction of HPßCD self-association and interaction with water. We then use the new set of parameters to characterize the behavior of HPßCD toward the different amino acids. We observe pronounced interactions of HPßCD with both polar and nonpolar moieties, with a special preference for the aromatic rings of tyrosine, phenylalanine, and tryptophan. Interestingly, our simulations further highlight a preferential orientation of HPßCD's hydrophobic cavity toward the backbone atoms of amino acids, which, coupled with a favorable interaction of HPßCD with the peptide backbone, suggest a propensity for HPßCD to denature proteins.

The Role of Cyclodextrins against Interface-Induced Denaturation in Pharmaceutical Formulations: A Molecular Dynamics Approach.

Protein-based pharmaceutical products are subject to a variety of environmental stressors, during both production and shelf-life. In order to preserve their structure, and, therefore, functionality, it is necessary to use excipients as stabilizing agents. Among the eligible stabilizers, cyclodextrins (CDs) have recently gained interest in the scientific community thanks to their properties. Here, a computational approach is proposed to clarify the role of ß-cyclodextrin (ßCD) and 2-hydroxypropyl-ß-cyclodextrin (HPßCD) against granulocyte colony-stimulating (GCSF) factor denaturation at the air-water and ice-water interfaces, and also in bulk water at 300 or 260 K. Both traditional molecular dynamics (MD) simulations and enhanced sampling techniques (metadynamics, MetaD) are used to shed light on the underlying molecular mechanisms. Bulk simulations revealed that CDs were preferentially included within the surface hydration layer of GCSF, and even included some peptide residues in their hydrophobic cavity. HPßCD was able to stabilize the protein against surface-induced denaturation in proximity of the air-water interface, while ßCD had a destabilizing effect. No remarkable conformational changes of GCSF, or noticeable effect of the CDs, were instead observed at the ice surface. GCSF seemed less stable at low temperature (260 K), which may be attributed to cold-denaturation effects. In this case, CDs did not significantly improve conformational stability. In general, the conformationally altered regions of GCSF seemed not to depend on the presence of excipients that only modulated the extent of destabilization with either a positive or a negative effect.